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I examine spectral properties of a dissipative chaotic quantum map with the help of a recently discovered semiclas- 
sical trace formula. I show that in the presence of a small amount of dissipation the traces of any finite power of the 
propagator of the reduced density matrix, and traces of its classical counterpart, the Frobenius-Perron operator, are 
identical in the limit of Ti — > 0. Numerically I find that even for finite fi the agreement can be very good. This holds in 
particular if the classical phase space contains a strange attractor, as long as one stays clear of bifurcations. Traces of 
the quantum propagator for iterations of the map agree well with the corresponding traces of the Frobenius-Perron 
operator if the classical dynamics is dominated by a strong point attractor. 



I. INTRODUCTION 
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The interplay between chaos, quantum mechanics and dissipation is rather complex and the subject of strong 
current research activities (|-^|. The definition of chaos in classical mechanics via exponentially fast spreading trajec- 
tories can not be applied to quantum mechanical systems, since the notion of a trajectory does not exist in quantum 
mechanics. On a quantum mechanical level chaos manifests itself in the statistical properties of the eigenenergies and 
eigenfunctions. In the case of Hamiltonian systems the eigenenergies and eigenfunctions obey the universal statistics 
of large random hermitian matrices restricted only by general symmetry requirements like invariance under time and 
spin reversal [0J|]. While no rigorous proof of this conjecture is known yet, overwhelming numerical and experimental 
evidence has been accumulated 

Dissipation has at least two very important effects. The classical dynamics is altered profoundly. It is no longer 
restricted to a shell of constant energy in phase space, but phase space volume shrinks if no external source com- 
pensates for the energy dissipated. In a chaotic system with external driving dissipation typically leads to a strange 
attractor in phase space, i.e. a multi-fractal structure that is invariant under the dynamics and which has a dimension 
strictly smaller than the dimension of the phase space. The second effect is of quantum mechanical nature and more 
subtle: Dissipation destroys very efficiently the quantum mechanical phase information. This typically happens on 



time scales much shorter than classical ones and even with very tiny amounts of dissipation |12-14|. Therefore the 
system behaves more classically, and one might expect to find classical manifestations of chaos again. It was indeed 
shown in a variety of examples that appropriate quantum mechanical counterparts of the classical phase space density 
(like Husimi functions or Wigner distributions) approach a smeared out version of the strange attractor M l^JI^] . At 
the same time one might ask whether spectral properties approach their classical counterparts as well. This paper 
1-^ ■ shows that for certain spectral properties the answer is "YES" , even though the structure of the spectrum can be 
Q \ very different in both cases. 



My analysis is based on a recently discovered trace formula for dissipative systems which, in the spirit of Gutzwiller's 
celebrated formula |l7^8|, expresses traces of the propagator of the density matrix in terms of classical periodic orbits 



[M . I show in section [II that to lowest order asymptotic expansion in h the traces agree with the traces of the classical 
Frobenius-Perron propagator of the phase space density |l!| . 
In the next section I briefly review basic properties of dissipative quantum maps, semiclassical theory and the trace 
formula. In section |iy] I apply the trace formula to a dissipative kicked top and compare with numerical results for 
finite h. The main results are summarized in section [v|. 



II. DISSIPATIVE QUANTUM MAPS 



A. General remarks 



Dissipation is introduced on a quantum mechanical level most rigorously by the so-called Hamiltonian embedding 
plf . The system of interest is considered as part of a larger system including the "environment" to which energy 
can be dissipated. The total system is assumed to be closed so that it is adequately described by a Schrodinger 
equation. The degrees of freedom of the "environment" remain unobserved. The system at interest is described by 
a density matrix in which the environmental degrees of freedom have been traced out, usually termed the reduced 
density matrix pit). 



Dissipative quantum maps are maps of the reduced density matrix from a time t to a time t + T: pit + T) = Ppit). 
They are analogous to non-dissipative quantum maps, where the state vector of the system is mapped with a unitary 
Floquet matrix F, \ip(t + T)) — F\ip(t)}. In the dissipative case P is not a unitary operator and therefore has 
eigenvalues inside the unit circle. 

Maps, dissipative or not, are a natural description of a time evolution if an external driving of the system is periodic 
in time with period T. They give a stroboscopic picture which suffices if the evolution during one period is of no 
interest. Systems that are periodically driven are capable of chaos even if they have only one degree of freedom. I 
will restrict myself in the following to such cases. 

The maps that I consider are particularly simple in the sense that the dissipation is well separated from a remaining 
purely unitary evolution where the latter by itself is capable of chaos. The unitary part will be described by a Floquet 
matrix F acting on the state vector, so that the unitary evolution takes the density matrix from pit) to p'{t) = Fp{t)F^ . 
After the unitary part a dissipative step follows which I will describe by a propagator D. It takes the density matrix 
from p'(t) to pit + T) = Dp' it). The total map therefore reads 

pit + T)=D{Fp{t)F^) = Ppit). (2.1) 



Such a separation into two parts is not purely academic. A most obvious realization of (2.1) is given when a Hamilto- 
nian Hit) leading to the unitary evolution and the coupling to the environment can be turned on and off alternatively. 
This should be realizable for instance with atoms flying through a series of cavities where in each cavity either the 
unitary evolution or the dissipation is realized. Another example might be a billiard, in which the particle only 



dissipates energy when hitting the walls. But even if the dissipation cannot be turned off, the map (2.1) may still be 
a good description. For instance, if the dissipation is weak and if the entire unitary evolution takes place during a 
very short time, dissipation may be negligible during that time. This is the case if the entire unitary evolution is due 
to a periodic kicking. The dissipation can then be considered as a relaxation process between two successive kicking 
events. Finally, a formal reason for such a separation can be given when the generators for the unitary evolution and 
the dissipation commute. 

These ideas might become clearer with a particular model. Let me therefore introduce as prime model system a 
dissipative kicked top. 



B. A dissipative kicked top 

The dynamical variables of a top p2| , pl| are the three components J x ,y,z of an angular momentum J. I will only 
consider dynamics (including the dissipative ones) which conserve the absolute value of J, J 2 = j(J + 1) = const. In 
the classical limit (formally attained by letting the quantum number j approach infinity) the surface of the unit sphere 
lim J -_ i . o0 (J/j) = 1 becomes the phase space, such that one confronts but a single degree of freedom. Convenient phase 
space coordinates are 

p = J Z J J — cos 9 — p and tp = q, (2-2) 

where the polar and azimuthal angles 9 and ip define the orientation of the angular momentum vector with respect 
to the J x , J v and J z axes. The parameter J is defined as J = j + 1/2 and allows for more convenient expressions of 
most semiclassical quantities. Due to the conservation of J 2 the Hilbert space decomposes into (2j + 1) dimensional 
subspaces. The quantum dynamics is confined to one of these according to the initial conditions. The semiclassical 
limit is characterized by large values of the quantum number j which can be integer or half integer. Since the classical 
phase space contains (2j + 1) states, Planck's constant may be thought of as represented by 1/J. 

Consider a unitary evolution generated by the Floquet matrix 

F = c - i TJ J *e- i P J y . (2.3) 

The corresponding classical motion first rotates the angular momentum by an angle (3 about the y-axis and then 
subjects it to a torsion about the z-axis. The latter may be considered as a non-linear rotation with a rotation angle 
given by the J z component of J. The dynamics is known to become strongly chaotic for sufficiently large values of 
k and /?, whereas either k = or [3 — lead to integrable motion |ll[] . For a physical realization of this dynamics it 
might be best to think of J as a Bloch vector describing the collective excitations of two-level atoms, as one is used 
to in quantum optics. The rotation can be brought about by a laser pulse of suitably chosen length and intensity, and 
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t he torsion by a cavity that is strongly off resonance from the atomic transition frequency [£3|. The Floquet matrix 
( |2.3|) has also been realized in experiments with magnetic crystallites with an easy plane of magnetization |24f| . 
Our model dissipation is defined in continuous time r by the Markovian master equation 



^P(r) = ^j([J-, P (r)J + ] + [J-p(r), J+]) = Ap(r) . 



(2.4) 



where the linear operator A is defined by this equation as generator of the dissipative motion. Equation (2.4) is 
well-known to describe certain superradiance experiments, where a large number of two-level atoms in a cavity of bad 
quality radiate collectively p5| , p6| . The angular momentum operator J is then again the Bloch vector of the collective 
excitation and the J+, </_ are raising and lowering operators, J± = J x ± U y . One easily verifies that ( |2.4| ) conserves 
the skewness in the J z basis (J z \m) = m\m)), i.e. matrix elements {m\\-j- p\m,2) with a given skewness Jrj — ni\ — ni2 
depend only on matrix elements with the same skewness. Eq.(2.4) is formally solved by p(r) = exp(Ar)p(0) for any 



initial density matrix p(0) and this defines the dissipative propagator 

D{t) = cxp(Ar) . (2.5) 

Explicit forms of D can be found in [p5|j2^j2^1 . The skewness rj only enters as a parameter in D. The classical limit 
gives the simple picture of the Bloch vector creeping towards the south pole 9 = n as an over-damped pendulum, 
according to the equations of motion 



^-6 = sin6>, -f-w = 0. 
(it or 



(2.6) 



Classically the azimuth tp is therefore conserved. Eq.(2.6) also shows that r is the time in units of the classical time 
scale. In the following it will be set equal to the time between two unitary steps. 



The Floquet matrix (|2.3|) is usually generated by the Hamiltonian 



H(t) 



2JT 



OO 

£ 
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S(t - nT) 



(2.7) 



it describes the evolution from immediately before a kick to immediately befo re the next kick. The g ene rator (p^j) for 
the dissipation does not commute with H(t). In order to obtain the map (2J) one should replace in ( |2.7| ) Ho = 2§T^ 
by 2 jTi J* anc ^ switch on H{t) only for a time T\ < T during each period T, whereas the dissipation acts during 
the rest of the time r = T — T\. Alternatively, when H(t) and A act permanently one may go to an interaction 
representation by p(t ) = exjp(—j-Hot)p(t) exp(j-HQt). In the J 2 representation this leads only to phase factors in the 
master equation (2A) which can be easily incorporated in D and which vanish moreover f or d ia gona l elem ents . Let 
us assume in the following that either has been done and use (2.1) with F and D given by (2.3), (|]J|), and (2.5) as a 



starting point with r as fixed parameter that measures the relaxation time between two unitary evolution and thus 
the dissipation strength. 



C. The trace formula 



In 1970 Gutzwiller published a trace formula for Hamiltonian flows that has become a center piece of subsequent 
studies of quantum chaotic systems [[l^,[l8| . A corresponding formula was obtained later for non-dissipative quantum 
maps by Tabor [^9). Assuming the existence of a corresponding classical map y = f (x) of phase on itself (x = (p', q') 
are the old, y = (p, q) the new phase space coordinates), both formulae express a spectral property of the quantum 
mechanical propagator as a sum over periodic orbits of f. Each periodic orbit contributes a weighted phase factor, 
where the weight depends on the stability matrix M of the orbit and the phase is basically given by the classical 
action S in units of h. Tabor's formula aims at traces of the Floquet matrix F, 

trF = Y, |2^teMjT72 ' < 2 - 8 > 

p.p. 1 1 

I have written the sum over periodic orbits as a sum over periodic points (p.p.) of the N times iterated map f N ; 
the integer v (the so-called Maslov index) counts the number of caustics along the orbit. All quantities have to be 
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evaluated on the periodic points. The squared modulus of trF N has in the unitary case an interpretation as (discrete 
time) form factor of spectral correlations. 



In [BO] a corresponding trace formula for dissipative quantum maps of the form (2.1) was derived. It is based on 
semiclassical approximations for both F and D. The semiclassical approximation of F has the general form of a van 
Vleck propagator |3(],|3l]]; a corresponding semiclassical approximation for D was obtained in |2q| . A WKB ansatz 
lead to a fictitious Hamiltonian system which depends on the skewness as a parameter. Its trajectories connect initial 
and final points specified by the arguments of D. Much as in the unitary case, an action R is accumulated along the 
trajectories; it has the usual generating properties of an action. Based only on the general van Vleck forms of F and 
D and the generating properties of the actions S and R we derived the trace formula 

E= *** 

trpjV =E | ~, " w N = l,2,.... (2.9) 

P . P . \txUU N M^ -tvM 

The sum is over all periodic points of the iV-times iterated dissipative classical map; the Ri are the actions of the 
fictitious Hamiltonian system for vanishing skewness accumulated during the ith dissipative step. The denominator 
contains the stability matrices Afl 1 ' for the it h di ssipative step and M for the entire map f N . The matrix with 



index i = N is a t the left of the product. Eq.(2i)) is a leading order asymptotic expansion in 1/J for propagators P 



of the type fl2.ip . The following restrictions apply: 

• The phase space is two dimensional. 

• The classical limit of the dissipative part of the map conserves one phase space coordinate (the azimuthal 
coordinate ip for the dissipation described by ( |2 . 4| ) ) . 

• The propagator D for the dissipative part conserves the skewness 77 of the density matrix in a suitably chosen 
basis and D has a single maximum as a function of 77 at 77 = 0. As indicated the dissipation (^4|) conserves the 
skewness in the J, basis. 



• The dissipation exceeds a certain minimum value. It is given by r ^ 1/J for the dissipation (2.4) and thus may 
become infinitesimally small in the classical limit J — - > 00. 



Eq.(2.9) shows that periodic orbits and classical quantities related to them still determine the spectral properties 



of the quantum system even in the presence of dissipation. The formula will now be studied in more detail. 



III. CONNECTION TO CLASSICAL TRACE FORMULA 



Remarkable about (2.9) is its simplicity. First of all, when propagating a density matrix, one would expect a double 
sum over period ic p oints. Indeed, in the dissipation free case one easily shows trP = jtrJ"! 2 , and trF is given by the 



Tabor formula (2.S) as a simple sum over periodic points |29j. Out of the double sum, only the "diagonal parts" 
survive. Decoherence induced through dissipation destroys the interference terms between different periodic points. 
For the diagonal terms the actions S and —S stemming from F and F' and the phases due to the Morse indices cancel. 
The square roots in the denominator combine to a power 1. Due to the cancellation of the phase factors the traces 



(2.9) are always real and positive. They fulfill herewith a general requirement for all propagators of density matrices 
that follows from conservation of positivity of the density matrix. On the other hand one may wonder whether the 
trace formula should not be an entirely classical formula, if all interference terms are destroyed. This is indeed what 
I am going to show now. 

The classical propagator of phase space density is given by P c /(y,x) = 6(y — f(x)). In the case where P c i{y, x) 
describes the map arising from the evolution during a finite time of an autonomous system, P c i is commonly called 
the Frobenius-Perron operator. For brevity I use the same name in the present dissipative situation. The trace of 
the Nth iteration of P c i is given by [fl9| 

^=EE^% (3-D 

p.o. r—1 1 v r ' 1 

= J2 r^J ,nn (3-2) 



|dct(l - M)\ ' 

p.p. 1 v ' 1 
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where the first sum in (3.1) is over all primi tive periodic orbits of length n p , r is their repetition number and M p the 
stability matrix of the primitive orbit. In (3.2), p.p. labels all periodic points belonging to a periodic orbit of total 
length N, including t he r epetitions, and M is the stability matrix for the entire orbit. 

The fact that M in (2.9) is a 2 x 2 matrix leads immediately to det(l — M) = 1 + detM — trM. Since the map 
is a periodic succession of unitary evolutions (with stability matrices M„^) and dissipative evolutions (with stability 



matrices M^ ), M is given by the product M = JX} = 
det m£ ] = 1 for alH = 1 . . . N and det M = J] I 



M®M®. 



conserves q which means that is diagonal, 



=N detM®. 



(i) 

The stability matrices Mi are all unitary so that 



The dissipative process for which (2.9) was derived 



M$> = 



1 

n ™« 



(3.3) 



The upper left element is 



M (i) 

i—N lvl d 



trM 



Sq(p',g') 

dq' '■ 



the lower right 



9p(p',q') 



But then detM 



(0 



and we find with 



= |1 + det M - trM | = | det(l - M)| exactly the denominator in (BJ). 



trn- 

The actions Ri are zero on the classical trajectories for the dissipative process (2.4), as one immediately sees by using 
their explicit form p8[ . Their vanishing can be retraced more generally to conservation of probability by t he master 
equation and therefore holds for other master equations of the same structure as well. To see this write (^J) in the J z 
basis and look at the part with vanishing skewness, i.e. the probabilities p m = (m\p\m). We obtain a set of equations 



dr 



Pm = {gm+lPrn + 1 - 9mPm) 



(3.4) 



where the specific form of the coefficients g m is of no further concern. Important is rather that the same function 
g m appears twice. This is sufficient and necessary for the conservation of probability, trp = Y^ m =-jP m = 1- Ori the 
other hand, had we coefficients f m and g m (i.e. p m = {g m +iPm+i ~ fmPm)) we would obtain the action R on the 
classical trajectory as JR — X)™=m( m (5'i) — m (/0) as one easuv verifies by writing down the exact Laplace image of 
D following the lines in p7| . Thus, the act ion is zero iff probability is conserved. But then the trace formula ft2.d[ ) is 
identical to the classical trace formula (3.1). 

This result proves that the traces of any finite power of the evolution operator of the quantum mechanical density 
matrix, are, in the limit of h~ — > 0, exactly given by the corresponding traces of the evolution operator of the classical 
phase space density, provided a small amount of dissipation is introduced. This is quite surprising since it is clear that 
even the basic structure of the two spectra can be very different: For all finite Hilbcrt space dimensions d = 2j + 1 
the quantum mechanical propagator P can be represented as a finite d 2 x d 2 matrix. Its spectrum is therefore always 
discrete, regardless of whether the corresponding classical map is chaotic or not. On the other hand, it is known that 
P c i has necessarily a continuous spectrum if the classical dynamics is mixing [|32"| . 

A formal reason why the spectra may differ in spite of the fact that the traces agree to lowest order in h is easily 
found. In order to construct the entire spectrum of P one needs d 2 = (2j + l) 2 traces. But already for traces of 
order j the next order corrections in the asymptotic expansion in 1/J that lead to (2.9) become comparable to the 
classical term; and for the highest traces needed (i.e. traces of order j 2 ) the next order in 1/J would be even more 

important than the classical term, so that one may not expect tvP 3 — trP^ for j — > oo. In other words, if we do not 
keep n fixed in the classical limit, trP" = trP c ™ may not hold for j — > oo and therefore the spectrum of P can be very 
different from that of P c ;. 



The asymptotic equality of trP" and trP^ strengthens substantially the quite envolved semiclassical derivation of 



(tL3> m 

should t 



It also sheds l ight on the question what happens if the dissipation does not conserve the coordinate q. We 
len not expect (2.9) to be valid but presumably replace it with the more general form (3.2). 



IV. COMPARISON WITH NUMERICAL RESULTS 



The question arises how good the agreement between quantum and classical traces is for finite J. To answer 
this question I have calculated numerically the exact quantum mech anical traces for our dissipative kicked top and 
compared them with the traces obtained from the trace formula (^^). These results will be presented now. 
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A. The first trace 



The quantum mechanical propagator P is most conveniently calculated in the J z basis, since the torsion part is then 
already diagonal. The rotation about the y-axis leads to a Wigner <i-function whose values are obtained numerically 
via a recursion relation as described in fl30| . The propagator for the dissipation is obtained by inverting numerically 
the exactly known Laplace image p5| , p7| . The total propagator P is a full, complex, non-hermitian, and non-unitary 
matrix of dimension (2j + l) 2 x (2j + l) 2 . Since for the first trace the knowledge of the diagonal matrix elements 
suffices I was able to calculate trP up to j = 80. Higher traces are most efficiently obtained via diagonalization which 
limited the numerics to j < 40. 

The effort for calculating the first classical trace is comparatively small. In all examples considered and even in the 
presence of a strange attractor, P c i had at most 4 fixed points that could easily be found numerically by a simple 
Newton-method in two dimensions. For each fixed point the stability matrix is found via the formulae in Appendix 
^ and so the trace is immediately obtained . 

In Fig.0 I show trP for different values of j as a function of r and compare with trP c /, eq.(2.9). The values for 



torsion strength and rotation angle, k = 4.0 and (3 = 2.0 were chosen such that the system is already rather chaotic 
in the dissipation free case at r = 0; a phase space portrait of many iterations of P c ; shows a large chaotic sea and 6 
relatively small stable islands. When t reaches a value of the order r ~ 0.5 a strange attractor appears which rapidly 
changes its form and dimension when r is increased. The attractor shrinks and is pushed more and more towards the 
south pole, as the angular momentum has more and more time to relax towards the ground state J z — —j between 
two kicks. At values of r of the order of r ~ 2.0 — 3.0 the attractor degenerates to a strong point attractor close to the 
south pole which absorbs even remote initial points in very few steps. At even stronger damping the point attractor 
reaches the south pole asymptotically. 

Figure |l| shows that - with the exception of very small damping - trP c ; reproduces trP perfectly well for all r, in 
spite of the strongly changing phase space structure. The agreement extends to smaller r with increasing j, as is 



to be expected from the condition of validity of the semiclassical approximation, r > 1/ J 28 1 . The analysis of the 
fixed points shows that at k = 4.0, (3 = 2.0 always two fixed points exist for r > 0.1. Their n component slowly 
decreases and the lower one converges towards the south pole with increasing r, where it finally coincides with the 
point attractor. 

Fig.|| shows the fixed point structure for a more complicated situation (k = 8.0, (3 = 2.0). The dissipation free 
dynamics at r = is entirely chaotic, no visible phase space structure is left. The above statements about the creation 
of a strange attractor (see Fig.^) and its degeneration to a point attractor when r is increased apply equally well. 

In Fig.[| I show the first trace as function of r for this situation. The classical trace diverges whenever a bifurcation is 
reached. Such a behavior is well known from the Gutzwiller formula in the unitary case; the reason for the divergence 
is easily identified as breakdown of the saddle point approximation in the semiclassical derivation of the trace formula. 
Whereas the quantum mechanical traces for small j (say j ~ 10) seem not to take notice of the bifurcations, they 
approximate the jumps and divergences better and better when j is increased. At j — 80 the agreement with the 
classical trace is already very good between the bifurcations. Remarkable, however, is the fact that there are some 
values of r close to the bifurcations, where all trP curves for different j in the entire j range examined cross. The 
trace seems to be independent of j at these points, but they nevertheless do not lie on the classical curve. One is 
reminded of a Gibbs phenomenon, but I do not have any explanation for it. 



B. Higher Traces 

Let us now examine higher traces tiP N for given values of k, (3, and r as a function of N. For large N all higher 
traces must converge exponentially to 1, independent of the system parameters. This is due to the fact that P has 
always one eigenvalue equal to 1. Its existence follows from elementary probability conservation [ jljj . The correspond- 
ing eigenmode is an invariant density matrix, its classical counterpart the (strange or point) attractor, the fixed points 
or any linear combinations thereof [ p2| . All other eigenvalues have an absolute value smaller than 1 since there is only 
dissipation and no amplification in the system. Their powers decay to zero as a function of N. 

I will focus on two limiting cases: The case where the basic phase space structure is a point attractor and the case 
where it is a well extended strange attractor. As explained above a point attractor can always be obtained by suffi- 
ciently strong damping. Consider the example k = 4.0, f3 — 2.0 and r = 4.0. Fig|| shows that indeed both quantum 
mechanical and classical result converge rapidly towards 1, and the agreement is very good even for j = 10. If one 
examines the convergence rate one finds that it is slightly j-dependent, but rapidly reaches the classical value. It 
should be noted that the calculation of trP c ^ is enormously simplified here by the fact that with increasing N no 
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additional periodic points arise. The dissipation is so strong that the system is integrable again. In the example 
given there are only two fixed points, one at (//, <ft) ~ (—0.3812219,-3.098751), a strong point repeller, and one at 
(/j,,4>) ~ (—0.9984018, —1.444154) a strong point attractor, and all periodic points of P£[ are just repetitions of these 
two points. 

The situation is quite different in the case of a strange attractor (Fig.^). The number of periodic points increases 
exponentially with N, as is typical for chaotic systems. This makes the classical calculation of higher traces exceed- 
ingly difficult. For k = 8.0, (3 = 2.0, and r = 1.0 I was able to calculate trPj^ reliably up to N — 5, where about 400 
periodic points have to be taken into account. The obtained numerical result for trPj^ can always be considered as 
lower bo und for the exact result for trP c ^ as long as one can exclude over-counting of fixed points since all terms in 
the sum (2J3) are positive. It is then clear that at N = 5 the quantum mechanical result at j — 40 is still more than 
a factor 3 away from trPj^, even though for N = 1 the agreement is very good. The convergence of trP^ to trPj^ as 
a function of j becomes obviously worse with increasing N. 



V. SUMMARY 



I have shown for certain dissipative quantum maps that the traces of (iterations of) the propagator of the quantum 
mechanical density matrix agrees to first order in an asymptotic expansion in Ti with the traces of the classical 
Frobenius-Perron propagator of the phase space density if a small amount of dissipation is present. This holds in 
spite of the fact that the corresponding spectra are very different. I have tested the theory numerically for finite 
values of h for a dissipative kicked top and have found good agreement in parameter regimes that ranged from very 
weak to strong dissipation. The phase space structure turned out to be important in the sense that higher quantum 
mechanical traces agree with very high precision with the classical ones if the phase space is dominated by a point 
attractor (strong dissipation) , whereas the precision is lost for higher traces in the case of an extended strange attractor 
(weak dissipation) . Sufficiently far away from bifurcations the lowest traces always agree very well with their classical 
counterpart. 
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tality of the ICTP Trieste, where part of this work was done. Numerical computations were partly performed at the 
John von Neumann-Institute for Computing in Julich. 



APPENDIX A: CLASSICAL MAPS AND THEIR STABILITY MATRICES 



I give here the classical maps for the three components rotation, torsion and dissipation as well as their stability 
matrices in phase space coordinates. All maps will be written in the notation (/i, (f>) — > (v, tp), i.e. /i and v stand for 
the initial and final momentum, (f> and ip for the initial and final (azimuthal) coordinate. The latter is defined in the 
interval from — n to n. The stability matrices will be arranged as 

/ m au \ 

M= [M fu • (Al) 



dfj, d^i 



a. Rotation by an angle (3 about y-axis 
The map reads 

v = [i cos (3 — yl — /i 2 sin j3 cos <fi (A2) 
■0 = ( arcsin(y // i^^sin0)6l(.T / ) + (A3) 

(sign(0)7r — arcsin(W sin (j))9{— x')J mod27r (A4) 

x' = y/l — /i 2 cos 4> cos (3 + [i sin (3 , (A5) 

where x' is the x component of the angular momentum after rotation, 9{x) the Heaviside theta-function, and sign(a;) 
denotes the sign function. 

The stability matrix connected with this map is 
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M r 



t/sin -0(^/l-^ 2 cos g+jj cos cf> sin g) p sin cog g , p cos sin j3 

\ ^/l-/i 2 (l-i/ 2 )cos-i/> ^(l-i/ 2 )(l-^ 2 )cosi/> " y/l-IJ. 2 



b. Torsion about z-axis 

Map and stability matrix are given by 



c. Dissipation 

The dissipation conserves the angle 0, and the stability matrix is diagonal: 



/i — tanh t 




(A6) 



v = n (A7) 
■0 = (</> + fyt)mod27r (A8) 

M t =( l k \). (A9) 



(A10) 
(All) 

M d = ( n l-(tanhr) 2 J • (A12) 
(1-ptanh-r) 2 / 

The total stability matrix for the succession rotation, torsion, dissipation is given by M — MdM t M r . 
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FIG. 1. Comparison of quantum mechanical traces (j = 10 (circles), j — 30 (squares), j — 50 (diamonds) and j = 80 
(triangles)) with classical trace (thick dashed line) for k = 4.0, P = 2.0 as function of r. There are no bifurcations for r > 0.1. 
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FIG. 3. Strange attractor for k = 8.0, /3 = 2.0, r = 1.0. The diamonds mark the position of the two fixed points. The 
borders tp — it and tp — —n have to be identified. 
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FIG. 4. Comparison of quantum mechanical traces with classical trace for k — 8.0, /3 = 2.0 as a function of r (same symbols 
as in Fig.nl). The classical trace diverges whenever a bifurcation is reached. 
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FIG. 5. Quantum mechanical and classical trace as a function of N for fc = 4.0, (3 = 2.0, r = 4.0 (same symbols as in 
Fig.^). The classical trace is shown as dashed line for better visibility, even though it is only defined for integer N. The inset 
shows that the exponential convergence to 1 also holds in the classical case. The classical dynamics is dominated by a single 
point-attractor/repeller pair. 
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FIG. 6. Quantum mechanical and classical trace as a function of N for k = 8.0, /3 = 2.0, r = 1.0 (same symbols as in Fig 
The corresponding phase space portrait is the strange attractor shown in Figpl 
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